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Abstract 

We outline a general strategy for the non-perturbative renormalisation of composite opera- 
tors in discretisations based on Neuberger fermions, via a matching to results obtained with 
Wilson-type fermions. As an application, we consider the renormalisation of the four-quark 
operators entering the AS = 1 and AS = 2 effective Hamiltonians. Our results are an 
essential ingredient for the determination of the low-energy constants governing non-leptonic 
kaon decays. 
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1 Introduction 



The renormalisation of four-fermion operators is an essential ingredient in lattice 
QCD computations of weak matrix elements. In this letter we will address the 
logarithmically divergent renormalisation of left-left four-quark operators, with an 
emphasis on the AS = 1 effective Hamiltonian governing non-leptonic kaon decays. 

The treatment of AS = 1 weak decays via an effective weak Hamiltonian with 
an active charm quark has been recently reviewed in pQ. After performing the 
operator product expansion and neglecting top quark effects, which are suppressed 
by three orders of magnitude relative to the contributions of up and charm quarks, 
the expression found for the AS = 1 effective weak Hamiltonian in the formal 
continuum theory is 

n w = J^(v us yv ud £{*f o? + *£G5} • (i-i) 

In the above expression g w = A^GpM^, kf 2 are Wilson coefficients, and the 
dimension-six operators QJ 2 have the form 

Qf = [(57 M P_u)(u7 M -P-<0 ± (s^P-d)(u^P-u)] - [u -> c] , (1.2) 
Qf = (ml - m 2 c ){m d (sP+d) + m s (sP_d)} , (1.3) 

where parentheses around quark bilinears indicate colour and spin traces and P± = 
|(1 ± 75). Although our procedure is completely general, we will from now on 
concentrate in the SU(4)l x SU(4)r symmetric limit, where all quark masses are de- 
generate In this limit the only contribution to decay amplitudes comes from 
matrix elements of the operators Qf~. Moreover, the operator renormalisation 
pattern is greatly simplified, as mixing with lower dimension operators is absent. 
We stress, however, that our results, being obtained in a mass-independent renor- 
malisation scheme, will renormalise properly subtracted operators also beyond the 
SU(4)l x SU(4)r symmetric limit. 

Our strategy to renormalise Qf is similar to the technique proposed in [2] for 
the computation of the renormalised chiral condensate. It involves matching bare 
correlation functions (or matrix elements) computed with Neuberger fermions to 
their renormalisation group invariant (RGI) counterparts. The latter are computed 
in the continuum limit with some variant of Wilson fermions, for which mature tech- 
niques for fully non-perturbative renormalisation exist. Our choice will be twisted 
mass QCD (tmQCD) with an 0(a) improved fermion action. 

Although we will concentrate specifically on the operators of the AS = 1 Hamil- 
tonian, both the proposed methodology and our results have a wider range of appli- 
cation. In particular, the renormalisation factors that we will obtain for Qf renor- 
malise also the four-fermion operator entering the AS = 2 effective Hamiltonian. In 
the present work, all computations are performed in the quenched approximation. 
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In the next section we will describe the strategy of the computation. In section 3 
we discuss the computation of RGI operators in the continuum limit, based on a 
twisted mass QCD (tmQCD) Wilson fermion regularisation. In section 4 we discuss 
our results with Neuberger fermions, and compute non-perturbative renormalisation 
factors. Section 5 deals with perturbative estimates of the same renormalisation 
factors. We present our conclusions in section 6. 



2 Strategy 

Let us consider a generic multiplicatively renormalisable operator. 1 The notation 
will follow closely that of pQ. We will be dealing only with mass-independent renor- 
malisation schemes. 

We start by recalling the definition of renormalisation group invariant (RGI) 
composite operators. The RGI insertion of a local operator Q into a continuum 
on-shell correlation function is given by 

Q(9R,m R: L) =c s (/i/A) lim Z s (g ,afi) Q(g ,m ,L) 

. (2.1) 
= lim Z (g ) Q(g ,m ,L) , 

where Z s is a renormalisation constant that renders the operator finite, A is the 
QCD scale, and go, mo (9R, m R) denote the bare (renormalised) gauge coupling and 
quark mass(es). The subscript "s" labels the renormalisation scheme. We have also 
indicated explicitly that correlation functions will be computed in a finite volume of 
spatial size L (eventually taking L — > oo). The RG-evolution function is given by 



c^/A) = [2&o3 2 (^)] 7 ° /(2M exp J - J'™ dg 



7(g) 7o_ 
P(g) b g 



(2.2) 



where we have used the perturbative expansions of the anomalous dimension of the 
operator 7 and the /3-function, viz. 

Pig) -g\bo + hg 2 + ...), 7 ( 5 ) 9 Z° g 2 ( l0 + iig 2 + ■■■)■ (2.3) 

It has to be stressed that the (scale-independent) RGI renormalisation factor Z(go) = 
Cg(/x/ A)Z s (<7o, a/j.) depends on the renormalisation scheme only via cutoff effects, 
since the RGI operator insertion Q is scheme-independent. On the other hand, Z is 
regularisation-dependent. We also stress that the running factor c s is a continuum 
quantity, and hence regularisation- independent. 

We now consider two different lattice regularisations, namely Wilson (denoted 
by "w") and Neuberger (or overlap) fermions (denoted by "ov"). Our aim is to 

1 The generalisation to operators which mix under renormalisation is straightforward. 
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construct RGI renormalisation factors for Neuberger fermions operators by matching 
renormalised quantities in both regularisations. The first step consists of using the 
first regularisation in order to compute the RGI operator Q at a reference physical 
point, parametrised here by (#r, ref , m^ref > L Ie f), viz. 

S(#R,ref > ™R,ref , ^ref) = Z w ' (g ) Q W {g , m , L re f) . (2.4) 

It is essential to note that any reference to the regularisation employed in the r.h.s. 
of eq. Q2.4J) has disappeared after the continuum limit has been taken. The second 
step consists of tuning a point (g'o,m' ) in the bare parameter space of the second 
regularisation, which corresponds to the same values of the renormalised parameters 
(#R ref > m R ref)- Assuming universality of the continuum limit, one then has 

QG?R,ref , i%ef, ^ref) = Z ov (g' ) Q ov (g' ,m' , L Iei ) + 0(a 2 ) , (2.5) 

where we have explicitly used the fact that correlation functions computed with 
Neuberger fermions exhibit scaling violations of at most 0(a 2 ). Once the bare 
quantity Q° v (g' , m' , L re f) has been computed, eq. Q2.5JI yields Z ov (g' ), provided 
that Q(gR re f , mR^ef, L re f) has been determined through eq. (|2.4|) . This procedure 
can be repeated at several bare couplings and masses (g' ,m' ), always corresponding 
to (<7r re f > m R,ref); m this wav the lattice spacing may be varied, while the physics 
(i.e. physical volume and renormalised coupling and masses) is kept fixed. Note that 
Eq. Q2.5|) is to be interpreted as a renormalisation condition that implicitly defines a 
mass-independent renormalisation scheme. Thus it ensures that the RGI renormal- 
isation factors computed in this way will correctly renormalise the operators at any 
value of the quark masses ttir. In particular, Z ov depends on quark masses only via 
cutoff effects (though it is not guaranteed a priori that such dependence is small). 
On the other hand, the renormalisation prescription l|2.5|) reproduces by construc- 
tion the RGI result at the reference point for any chosen set of bare parameters. 
Thus the procedure is only useful if the targeted physical regime, characterised by 
{9R,m R ,L), is well away from (#R jref , m Rirc f , L ie{ ). 

The present work provides an application of this strategy. The ultimate aim, 
which is achieved in ref. [2] , is the computation of the effective low-energy couplings 
governing non-leptonic kaon decays, following the strategy described in pQ. This 
involves the computation, carried out using Neuberger fermions, of the chiral limit 
values of the ratios of correlation functions 

N±( s = ([Mx)] du Qt(o)[My)}us) , , 

1 1 ([Mx)] du [M0)} ud )([J (0)} su [J (y)} us ) ' 1 • ] 

where is the left-handed current 

[J^{x)] a (3 = ip a ^P-ipp , (2.7) 
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and a, (3 are flavour indices. An essential ingredient of the procedure are the renor- 
malisation factors [Z 1 /Zj\, the computation of which is the goal of the present 
work. 

In order to compute non-perturbatively the renormalisation factors [Z^/Z%] for 
Neuberger fermions, we will employ the ratios of QCD matrix elements 

{ir + \Qt\K+) 

U± = <^|J |0><0|Jo|K+> (2 ' 8) 

computed in large volumes and at a value of the reference quark mass m^ re f 
corresponding to mps = m n = = mF^ ys = 495 MeV. Note that, in the 
SU(4)l x SU(4)r symmetric limit, the ratio TZ+ will be equal, up to a trivial fac- 
tor, to the Kaon bag parameter Bk- The RGI ratios 1Z± will be computed, as in 
Eq. (|2.4|) . using a Wilson fermion regularisation. To this purpose, we will compute 
the bare quantities 1Z± at several values of the bare coupling go, and use the RGI 
renormalisation factors computed in the same go range in [1], using a Schrodinger 
Functional (SF) framework. We will then apply Eq. 112.5(1 to match the RGI ratios 
TZf to the ratios of bare matrix elements computed with Neuberger fermions. Since 
the matching reference regime of large volumes and meson masses of the order of 
m^ ys is well away from the target one in which [Z^/Z^] are to be used, the con- 
struction is indeed non-trivial. In this way we have exploited the fact that Wilson 
fermions are well suited for simulations in the strange-quark regime, while they be- 
come problematic close to the chiral limit, where Neuberger fermions are clearly 
advantageous. 

For the sake of consistency, we will also perform a direct determination of the 
ratio Z^~' ov /Z^' ov , computed from a matching involving the ratio of matrix elements 

K + (*+\Qt\K+) 

K- {n+\Q^\K+) ' y ' 1 

This specific ratio is of particular interest, as it enters directly the study of the 
AI = 1/2 enhancement rule. 

Some comments are in order. The renormalisation factors obtained via the 
procedure just described do not lead, obviously, to independent renormalised values 
for the ratios of matrix elements in Eq. (|2.8J) computed at the reference point m^ ys , 
as a tautology would result. As explained above, the renormalisation factors [Z^/Z]] 
will rather be used to renormalise quantities effectively computed in the chiral limit. 
They can also be used e.g. to renormalise ratios of correlation functions computed 
in the e-regime. The fact that the matching involves correlation functions computed 
with both periodic and SF boundary conditions does not, on the other hand, give 
rise to any subtlety, as only hadronic matrix elements computed in a large volume 
are involved. 
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Unlike the above strategy, adopted in this work, the ideal matching procedure 
should not involve the tradeoff of a long-distance matrix element of physical rele- 
vance. It is e.g. possible to match a different matrix element of the same operator, 
or to take a reference point for the matching which is well away from all the target 
physical regimes of interest. On the other hand, the particular strategy adopted 
here has the advantage that it allows to use the numerical results obtained in the 
context of |5;. 

3 Wilson-tmQCD computation of RGI operators 

We will now discuss the computation of the RGI operators Qf (i.e. the l.h.s. of 
Eq. (|2.5jl ). using the tmQCD formalism with Wilson quarks (Hj- 

We start by recalling that, with Wilson fermions, the renormalisation of Q 1 is 
more subtle than in chirally symmetric regularisations (see [Zj for a detailed discus- 
sion). Customarily, both operators are split into parity-even and parity-odd parts 
as 

Qf = Qvv+aa ~~ 2va+av ' 

in standard notation. In the three-point correlation functions considered below, par- 
ity conservation in QCD ensures that the only contribution comes from the parity- 
even part Qyv+AA- With ordinary Wilson fermions, as a consequence of the break- 
ing of chiral symmetry, the renormalisation of Qyv+AA requires the subtraction of 
four finite counterterms involving all the remaining Lorentz-invariant, parity-even 
four-quark operators with the same flavour structure. 2 On the other hand, it is pos- 
sible to construct twisted mass QCD (tmQCD) Wilson regularisations in which the 
counterterms are absent, and the operator renormalises multiplicatively. The basic 
property that has to be satisfied is that the chiral rotation of the quark fields that 
generates the twisted mass term maps Qyv+AA on t° QyA+AV ^ n mass- independent 
renormalisation schemes, the latter is protected from mixing with other four-fermion 
operators by CPS symmetry. Examples of such regularisations have been discussed 
in 0. 

Here we will adopt, however, a different approach, which will allow us to use 
the numerical results obtained in j^j. To that purpose we restrict ourselves to the 
quenched approximation, and use the formalism for valence quark flavours advocated 
in [0]. We consider a theory with six valence flavours, that we will label ip = 
(u, d, s, c, v! , c') T . We will employ two tmQCD regularisations, characterised by the 

2 Mixing with operators of lower dimension is always absent in the SU(4)l x SU(4)r limit, as all 
mixing coefficients are proportional to the mass difference (m c — m u ). 
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choice of twist angle a in the definition of the fermion action: 



StaQCD = « 4 E fa) { fl w,sw + mW + * 7 5M (a) } (x, y) ^(y) , (3.2) 

where -D w ,sw is the Wilson operator with a Sheikholeslami-Wohlert term, m( Q ) and 
fj,( a ^ are diagonal mass matrices, and the label a refers to the twist angle enter- 
ing chiral rotations. The bare mass parameters are tuned so that, up to 0(a 2 ) 
corrections, the renormalised mass matrices have the form 

m£ /2) = M R diag(0, 0,1, 0,1,1), /i£ /2) = Mr diag(l, 1, 0, 1, 0, 0) , (3.3) 
m £/4) = Mr d . ag(l] ^ ^ ^ ^ l} ^ /4) = Mr ^ _^ ^ _^ _^ ^ ^ 

where Mr is the physical renormalised quark mass. The mass tuning procedure is 
identical to the one described in We then introduce the operators 

Q± =[{s llx P_u){u'^P-d) ± (s^P„d)(u' lfl P„u)] - [u - c, u' - c'} . (3.5) 

Using the standard relation between QCD and tmQCD in the continuum limit, one 
finds, for the two regularisations specified above, that the following equalities hold 
between operator insertions in RGI correlation functions 

[QvV+Aa]qCD = i [QvA+AV"] tmQCD , (3.6) 
^a[^]sm,QCD = {Z A [Afj] SUj tmQCD ~ ^v[^/i]sw,tmQCD} , (3.7) 

Z A [Ap]du',QCD = \Z A [Aft] du' jtmQCD — iZy [V/j] du' , tmQCD } , (3.8) 

where [A^] a p = tpalfils^p and [V^] a(} = tpaln^p- The current normalisations 
Z A ,Z V , as well as the O(a) improvement coefficient c A needed to construct an O(a) 
improved axial current, are set to the values computed by the Alpha Collabora- 
tion dnnm. 

Correlation functions are computed within a Schrodinger functional (SF) frame- 
work, with quark and gluon fields obeying periodic boundary conditions in space 
(with period L) and (homogeneous) Dirichlet boundary conditions in time at the 
hypersurfaces xq = and xq = T. Ratios of correlation functions, from which the 
ratios of matrix elements 1Z± and 1Z+ /1Z— can be extracted, are defined in complete 
analogy to the ones specified for the extraction of Bk in Our run parameters, 
too, are the same as in and are listed in Tables 1 and 5 of that work, save for 
one important exception, concerning the dataset at (3 = 6.1. After completion of 
ref. [2j , its authors carried out a proper determination of K cr ((3 = 6.1), based on the 
method of ^U]. They found that this estimate of k ct disagrees by several standard 
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deviations from the value cited in the literature ^2], which had been used in 
for the tuning of the bare mass parameters. This discrepancy hence necessitated a 
new determination of the bare mass parameters in order to satisfy the constraints 
imposed by the prescribed values of the twist angles. Consequently, the runs at 
= 6.1 had to be repeated after completion of 5 . Full details will be provided 
separately [TI]. In the present work we merely quote the corrected value of 1Z+ in 
the continuum limit. 

The RGI ratios are obtained upon multiplication by the RGI renormalisation 
factors Zy A+Av . The latter have been computed non-perturbatively in [3], using 
the standard SF finite-size scaling analysis, for a range of inverse couplings 6.0 < 
(3 < 6.5. Out of the various SF renormalisation schemes considered in ^ we have 
chosen to employ scheme 1 for the renormalisation of QyA+AV an< ^ scheme 8 for 
Qva+AV' t ne reasons are explained in [3] and [Ej- The appropriate error analysis 
has been extensively discussed (for 1Z+) in The continuum limit is then obtained 
by performing a combined extrapolation of the results coming from both tmQCD 
regularisations. The extrapolation is linear in the lattice spacing a, since the four- 
fermion operator is not 0(a) improved, and hence the leading lattice artifacts in 1Z± 
are expected to be 0(a). Furthermore, as discussed in 0, at the lowest values of 
(3 the 0(a) ambiguity in the determination of the improvement coefficient c A has a 
significant impact on cutoff effects. Under these premises, our most stable continuum 
limit extrapolation for 1Z + is obtained by discarding the /3 = 6.0,6.1 data points, 
while for 1Z- and 1Z + /1Z- only = 6.0 is discarded. The final results, illustrated in 
Fig. d are 

TZ + = 0.885(86) , = 0.849(82) , (3.9) 

= 0.875(80) . (3.10) 

We stress that the volume dependence of these results is well within the quoted 
uncertainties (see [S] for details). 

Eqs. ()3.9l3.10j) are the main result of the present work. In the next section 
we will use them to determine the renormalisation factors needed with Neuberger 
fermions. It must be stressed at this point that the continuum limit extrapolation is 
rather long and, in the case of 7£_, strongly driven by the (3 = 6.45 datum. A better 
control of the continuum limit extrapolations could be achieved e.g. by removing 
O(a) effects as suggested in ref This is beyond the scope of the present letter. 



4 Renormalisation constants for Neuberger fermions 

Having constructed the RGI ratios of matrix elements in eq. (|2.8[) . we now insert 
them in Eq. ()2.5j) and solve for the Neuberger fermions renormalisation constants 
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p 


am tquips 








5.8485 


0.060 1.259(10) 


0.772(30) 


1.514(73) 


0.511(28) 



Table 1: Results with Neuberger fermions, obtained on a 16 3 x 32 lattice from 197 
configurations and using low-mode averaging with 20 low modes of the Neuberger- 
Dirac operator. The physical spatial extent of the lattice is L/(2tq) = 1.98. 



In order to regularise the theory using Neuberger fermions, we start by intro- 
ducing the Neuberger-Dirac operator 15 

D = ^{l-A(AlA)-V 2 } , A = l + s-aD w , (4.1) 

where D w is the massless Wilson-Dirac operator, a denotes the lattice spacing, and s 
is a free parameter in the range \ s\ < 1. By setting a = a/ (1 + s) it is straightforward 
to check that D satisfies the Ginsparg- Wilson relation 

75 D + D l5 = aD l5 D . (4.2) 

Composite operators which have proper chiral transformation properties in the reg- 
ularised theory are obtained by performing the substitution 

^-(l-aaD)V, 

(4.3) 

The operators in the discretised theory share the same transformation properties 
under chiral symmetry as their counterparts in the continuum (see [I] and references 
therein). 

Bare values for the ratios of matrix elements 1Z±, are extracted from the ratios 
of correlation functions of eq. (|2.6j) . The details of the computation, performed at 
a fixed value of (3 with periodic boundary conditions in all Euclidean spacetime 
directions, are reported in [T| 1161 13]. The simulation parameters and our results for 
1Z± are provided in Tabled Since our pseudoscalar mass is compatible within errors 
with the kaon mass rQin^ ys = 1.2544, there is no need to consider other values of 
the quark mass to inter/extrapolate the kaon point. 3 Again, finite volume effects 
are expected to lie within the quoted uncertainties. 

Finally, by combining the continuum limit results of Eqs. I|3.9l3.1()j) with the 
bare Neuberger fermions results of Tabled we derive non-perturbative estimates of 

3 The value of the reference scale ro is set to ro = 0.5 fm, and we take the ratio ro/a from |17) . 
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bare P.T. 


MFI P.T. 


non-perturbative 


zt/z\ 


1.242 
0.657 
0.525 


1.193 
0.705 
0.582 


1.15(12) 
0.561(61) 
0.584(62) 



Table 2: Perturbative and non-perturbative estimates for Neuberger fermions RGI 
renormalisation factors at = 5.8485. 



the RGI renormalisation factors 



/3=5.8485 



zi 



/3=5.8485 



K+/K- 



(4.4) 



The results are collected in the last column of Table El together with the correspond- 
ing perturbative estimates, which will be discussed in the next section. 



5 Perturbative estimates of renormalisation factors 

In this section we will determine the RGI renormalisation factors of interest in 
perturbation theory. This provides a handle on the systematics related to their 
non-perturbative determination. 

The anomalous dimensions ^ of the operators Qf are known at two loops 
for several schemes. For discretisations based on the Neuberger-Dirac operator, 
the renormalisation factors Z s (go,a/i) have been computed for s = RI/MOM in 
perturbation theory at one loop in [!§]. The ratios of renormalisation constants we 
are interested in, computed with Neuberger fermions and in the RI/MOM scheme, 
can be written as 



Z m (g ,a») + 3) _|^_ ( 21n(4H _ l (Bg _ By) \ + O(5o 4 } > 



Zligo) v y 167r 2 1 w ' 3 

Z f^\ = 1 + 4- 2 {12 ln(4/,a) - 2(5s - By)} + O( 5o 4 ) 



(5.1) 



It is also possible to perform the expansion using "mean-field improvement" (MFI) 
|19j . which aims at improving the convergence of the perturbative series. At the 
level of the ratios in Eq. ()5.1|) . it is easy to check that the implementation of MFI 
simply amounts to replacing the bare coupling g$ by a "continuum-like" coupling 
q 2 , which we set to be gf-^. 

The coefficients B§ and By in Eq. (|5,1[) are listed in Table 1 of 18]. In order to 
obtain the corresponding RGI renormalisation factors, it is enough to multiply the 
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above by the suitable perturbative running factors c RI (/i/A). In our simulations we 
use = 6/gl = 5.8485. For fj, = 2 GeV and A = 238 MeV the NLO perturbative 
values for the coefficients c^j are c^j(/i/A) = 0.6259 and Cpj(/i/A) = 1.2735. Putting 
this together with Eq. ()5.1|) . we obtain for the RGI renormalisation factors the 
values quoted in the first two columns of Table |2j It is worth mentioning that the 
differences between perturbative results evaluated in "bare" and MFI perturbation 
theory are relatively small. This is presumably a consequence of having considered 
ratios of operators, in which contributions of the self-energy type cancel, and is 
in stark contrast to the situation encountered in simple quark bilinears, where the 
deviations between perturbative and non-perturbative estimates amount to about 
30% at similar values of the bare coupling [2J 1211 122j . 

This analysis implies, furthermore, that it is unlikely that our non-perturbative 
results are affected by large cutoff effects, e.g. those proportional to powers of the 
quark mass. 

6 Conclusions 

In this work we have laid out a general strategy for the non-perturbative renormal- 
isation of operators with Neuberger fermions, via a matching to results obtained 
with Wilson-type regularisations. As an application, we have dealt with the overall 
logarithmic renormalisation of the operators entering the AS* = 1 effective Hamilto- 
nian with an active charm quark, for which we have computed RGI renormalisation 
factors in the quenched approximation. An immediate application of our results 
appears in the determination of the effective couplings governing kaon decays in the 
low-energy description of the theory [Sj , in the spirit of . 
There are a few caveats in this approach: 

• From the technical point of view, we believe that our tmQCD results for the 
RGI Q ± constitute a significant advance with respect to previous computa- 
tions, in that they have been achieved with two Wilson-type regularisations, 
non-perturbative renormalisation and RG running, at several bare couplings 
etc. In spite of this, the fact that continuum limit extrapolations are rather 
long renders the absence of 0(a) improvement an important drawback in 
our effort to obtain stable continuum limit results. As far as our Neuberger 
fermions computations are concerned, we point out that, at present, we have 
results only at one bare coupling. Furthermore, exploring the dependence of 
renormalisation factors on the choice of reference point would be useful to 
quantify the impact of 0((am) 2 ) cutoff effects. 

• From the conceptual point of view, the specific matching procedure adopted 
here is based on fixing the matrix element (irlQ^lK) (at m w = txik = rn v ^ ys ) 
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to the value predicted with tmQCD Wilson-type fermions. Having used this 
"physical" predictions as renormalisation conditions (for Q + it is the value 
of the kaon mixing parameter Bk) implies that our measurements of jZ\ 
cannot be used for the independent renormalisation of Bk or of K — > it matrix 
elements with Neuberger fermions. On the other hand, our renormalisation 
constants are perfectly suitable to renormalise K — * tttt matrix elements com- 
puted in infinite volume and for particle masses in the physical range, or for 
ratios of correlation functions computed in the e-regime of QCD. 

The ideal approach to the renormalisation problem in hand would involve a 
working formulation of the Schrodinger functional for Neuberger fermions. An im- 
portant recent step in that direction is the proposal of ref. |23j . 

The tmQCD data reported in this work were obtained within an ALPHA Col- 
laboration project; P.D., F.P., CP. and A.V. wish to thank J. Heitger and S. Sint 
for discussions. Our calculations were performed on the APEMille installation of 
DESY-Zeuthen and on PC clusters at DESY-Hamburg, CILEA and the University 
of Valencia, as well as on the IBM Regatta at FZ Jiilich and on the IBM MareNos- 
trum at the Barcelona Supercomputing Center. We thank all these institutions and 
the University of Milano-Bicocca (in particular C. Destri and F. Rapuano) for their 
support. P.H. acknowledges partial support by CICYT (grants FPA2004-00996 and 
FPA2005-01678) and Generalitat Valenciana (GV05-164). F.P. acknowledges the 
Alexander-von-Humboldt Stiftung for financial support. 
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Figure 1: Continuum limit extrapolations of 7^. ± and 
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